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We report the first evidence for the decay E + — » pfi + n~ from data taken by the HyperCP (E871) 
experiment at Fermilab. Based on three observed events, the branching ratio is £>(E + — > pfi + fi~) = 
[8.6^5'4(stat) ± 5.5(syst)] x 10~ 8 . The narrow range of dimuon masses may indicate that the decay 



proceeds via a neutral intermediate state, E + 
MeV/c 2 and branching ratio £(£+ -> pP°, P° 



pP°,P° 



with a P° mass of 214.3 ± 0.5 



) = [3.1Ii;|(stat) ± 1.5(syst)] x 10" 
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PACS numbers: 13.30.Ce, 14.20.Jn, 14.80.Mz 

In the standard model (SM), the decay S + — > 
pl + l~ (£pu,l = e,/i) can be described as proceeding 
through a flavor-changing neutral-current (FCNC) inter- 
action and by internal conversion, as shown in Fig.^a)- 
(c). Bcrgstrom et al. Q argue that in the SM the FCNC 
contribution for the decay E^~ ;i is not dominant. The de- 
cay is of interest since it also allows a direct search 
for a new scalar or vector particle, which could contribute 
an s — ► d transition at the tree level |2| (Fig. IHd)). 

Current literature reports only an upper limit B(T, + — > 
pe + e~) < 7 x 10~ 6 3]. The decay rate for the process 
S + — > pl + l~ was studied in 0, 0] in the context of the 
SM. Using the measured partial decay width and the a 7 
parameter for S + — » pj decay, Bergstrom et al. Q] find: 
r(E+ -> pe+e-)/r(S+ -> pi) > 7.2 x and 1/120 > 



T(E+ -» w + M -)/r(s+ 



pe + e 



) > 1/1210. A large 



violation of these limits would signal new physics. 

In this Letter we report the first evidence for the 
E + — > decay, a measurement of the branching 

ratio for this decay, and possible evidence for a new state 
P° observed via E+ -> pP°, P° -> /i+ fr . 

The HyperCP experiment was located in the Meson 
Center beamline at Fermilab. The spectrometer, shown 
in Fig. |21 is described in detail elsewhere Charged sec- 
ondary beams, with mean momenta of about 160 GeV/c, 
were produced by 800 GeV/c protons incident on copper 
targets and momentum selected by a curved collimator 



situated in a dipole magnet (hyperon magnet) . The sign 
of the charged secondary beam was periodically changed 
by reversing the field of the hyperon magnet. We an- 
alyzed 2.14 x 10 9 triggers from the positive-secondary- 
beam data set and 0.37 x 10 9 from the negative. 

The signature of the decay is two unlike-sign 

muon tracks and a proton track originating from a com- 
mon vertex. The transverse-momentum kick of the an- 
alyzing magnets was such that muons from this decay 
were always deflected to opposite sides and outside of 
the intense secondary beam at the rear of the spectrom- 
eter. Hence the signal trigger used to select candidate 
events required hits in the same-sign and opposite-sign 
hodoscopes (Left-Right trigger), in coincidence with hits 
in the vertical and horizontal hodoscopes in each of the 
muon stations situated on either side of the secondary 
beam. 
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FIG. 1: Feynman diagrams for decays in the SM (a)-(c) 
and via new physics (d). The SM processes are referred to as 
FCNC (a) and internal conversion (b)-(c). 
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FIG. 2: Plan view of the HyperCP spectrometer. 



The basic selection cuts discussed below were applied 
to the Si",,,, event candidates. Each event was required 
to have a track in each muon station and a higher- 
momentum track (the proton candidate) with the same 
charge sign as the secondary beam. The total momen- 
tum of the three tracks had to be within the range 120 
to 240 GeV/c. The S + trajectory had to extrapolate to 
within 3.5 mm (« 3 a) of the center of the target, where 
the nominal target position was determined using well- 
reconstructed K + — > 7r + 7r + 7r~ decays and the extrapola- 
tion resolution was determined from a Monte Carlo (MC) 
simulation of „ decays. The decay vertex of the three 
tracks was calculated by the method of distance of closest 
approach (DC A), and the z-coordinate of the vertex z v 
was required to be more than 68 cm downstream of the 
entrance of the vacuum decay region and more than 32 
cm upstream of its exit. The average distance between 
pairs of tracks in the x— y plane at z v was required to be 
less than 0.25 mm. The hits in the multiwire propor- 
tional chambers upstream of the analyzing magnets were 
refit with a constraint that they share a common vertex. 
The resulting x 2 /ndf was required to be less than 1.5. 
After imposition of the above cuts, three candidate 
„ events were observed in the positive-secondary- 
beam data with masses (M pflfl ) within la (1 MeV/c 2 ) 
of the S + mass, as shown in Fig.[!2a) and (b). No other 

mass. 



events were found within m 20 MeV/c 2 of the S H 

The signal events were verified by applying two addi- 
tional cuts that removed almost all of the higher-mass 
background events without affecting the three signal 
events. The first cut required that the ratio of the pro- 
ton momentum to the summed three-track momentum 
be larger than 0.68, as MC simulations showed that this 
cut preserved 100% of the signal while removing most 
of the K + decay background. The second cut removed 
events whose mass was within ±10 MeV/c 2 (3tr) of the 
K + mass when calculated using the ir + /x + pT decay hy- 
pothesis to reject K + — > 7r + ^ + /i~ decays. The resulting 
pfi + fj," invariant mass distribution is shown in Fig[3Jc). 

Possible background sources were extensively studied, 
using both MC and data. Typically 100 times as many 




FIG. 3: M PMM distribution for the positive-secondary-beam 
data, (a) after the standard cuts, (b) within ±15 MeV/c 2 of 
the E + mass, and (c) after the additional cuts. The arrow 
represents the E + mass. 



MC events as would be expected for each potential back- 
ground source were generated. 

Other positively-charged hyperon decays, S + — > 
A-7T+ — ► pn + ir + and fl + — > AK + — > pn + K + , were not 
significant backgrounds as decays in flight of the daugh- 
ter 7r + or K + would produce muons of only one po- 
larity, and the probability of misidentifying the sign of 
the charge of any of the decay daughters was negligibly 
small. Contributions from charged kaon decays such as 



K + — > 7T + 7r + 7T , K + — > 7T + 7T fJ, + Va, K + — > TT + /i + fl 



and K + — > [V fi + ^Tv^ were estimated using MC simu- 
lations allowing in-flight decay of pions, and found to be 
negligible. Muon pair production by photon conversion 
in material inside the vacuum decay region was stud- 
ied using MC simulation for K + — > 7r + 7r°,7r° — ► 77, 



K+ 



7T + 77, 



S + — > j»7, and S + — > p7r°,7r° 



77 de- 



cays. Such sources of background were also found to be 
negligible. 

In addition, we used real data to investigate possi- 
ble backgrounds otherwise missed, including the unlikc- 
sign dimuon sample from the negative-secondary-beam 
data, as well as a sample of events (single-muon sam- 
ple) for which only one muon track was required in ei- 
ther the left or right muon station. For the single-muon 
sample both the positive- and negative-secondary-beam 
data were included, and the non-muon track was re- 
quired to be within the fiducial volume of the appropriate 
muon station. For the positive-secondary-beam data, the 
single-muon sample was one order of magnitude larger 
than the unlike-sign dimuon sample. These background 
studies showed that after cuts none of these sources con- 
tributed in the pp + fi~~ invariant-mass region below 1200 
MeV/c 2 . Finally, we relaxed each cut to increase the 
background level shown in Fig. OJc) by an order of mag- 
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nitude. However, there still were no background events 
within 8 MeV/c 2 of the £+ mass. 

Figure a) compares the dimuon mass distribution of 
the three signal candidates with that expected in the 
SM with the form-factors described below. The recon- 
structed dimuon masses for the three candidates, 214.7, 
214.3, and 213.7 MeV/c 2 , all lie within the expected 
dimuon mass resolution of rj 0.5 MeV/c 2 . The dimuon 
mass distribution for Sp„„ decays is expected to be broad 
unless the form factor has a pole in the kinematically al- 
lowed range of dimuon mass. 

The expected SM distribution was used to estimate the 
probability that the dimuon masses of the three signal 
candidates be within 1 MeV/c 2 of each other anywhere 
within the kinematically allowed range. The probabil- 
ity is 0.8% for the form-factor decay model and 0.7% for 
the uniform phase-space decay model. The unexpectedly 
narrow dimuon mass distribution suggests a two-body 
decay, £+ -> pP°, P° -> p+ pT wh ere P° is an 

unknown particle with mass 214.3 ± 0.5 MeV/c 2 . The 
dimuon mass distribution for the three signal candidates 
is compared with MC Sip decays in Fig. IUV), and 
good agreement is found. Distributions of hit positions 
and momenta of the proton, [i + , and fj,~ of the three can- 
didate events were compared with MC distributions, and 
were found to be consistent with both decay hypotheses. 

To extract the branching ratio, the E + — » 

pir°, 7T° — > e+e~7 (S+ ee7 ) decay was used as the normal- 
ization mode, where the 7 was not detected. (HyperCP 
had no 7 detectors.) The trigger for the £p ee7 events was 
the Left-Right trigger prescaled by 100. The proton and 
two unlike-sign electrons were required to come from a 
single vertex, as were the three tracks of the signal mode. 

The proton was selected to be the positively-charged 
track with the greatest momentum, and the event was 
discarded if the proton candidate did not have at least 
66% of the total three-track momentum, as determined 
by a MC simulation of £+ ee7 decays. The reconstructed 
mass for the 37r hypothesis was required to be outside 
±10 MeV/c 2 of the K+ mass. The cuts on x 2 /ndf, DCA, 
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FIG. 4: Real (points) and MC (histogram) dimuon mass dis- 
tributions for (a) MC events (arbitrary normalization) 
with a form-factor decay (solid histogram) and uniform phase- 
space decay (dashed histogram) model, and (b) E+p MC 
events normalized to match the data. 



and the total momentum were the same as for the signal 
mode. However, the decay vertex had to be more than 
168 cm downstream of the entrance of the vacuum decay 
region and more than 32 cm upstream of its exit. Since 
the 7 momentum was not measured, the x and y positions 
of the E + trajectory at the target were determined using 
only the three charged tracks, and those positions had to 
be consistent with that expected from a MC simulation 
of £+ e7 decays. 

To significantly reduce contamination from photon- 
conversion events, the dielectron mass was required to 
be between 50 and 100 MeV/c 2 . After application of the 
above selection criteria, a total of 211 events remained, 
as shown in Fig. [SJ We performed a binned maximum- 
likelihood fit for the mass distributions for data and three 
MC samples: £p ee7 decays, K + — » 7r + 7r°, 7r° — * e + e~7 
(K+ ee7 ) decays, and uniform background. From the 
fit, the number of observed S+ ee7 decays was = 
189.7 ± 27.4 events, where the uncertainty is statistical. 
To extract the total number of normalization events, val- 
ues of (51.57±0.30)% and (1.198±0.032)% were used re- 
spectively for the E + — > pir° and ir° — > e + e~7 branching 
ratios 

The kinematic parameters for S + production at the 
target were tuned to match the data and MC £p~ ee7 mo- 
mentum distributions. The MC E+, decays were gen- 
erated using the decay model in Ref. |7j for tt° — ■> e + e~7 
(7Tg e7 ) decays, and the ir° electromagnetic form- factor pa- 
rameter a = 0.032 ± 0.004 was taken from Ref. Q. After 
tuning of the parameters, comparisons of the distribu- 
tions of the MC events with the data for S+ ee7 decays, 
the decay vertex positions, momentum spectra, recon- 
structed mass, hit positions of each charged particle, etc., 
showed good agreement. 

In the simulation of the £p MM decays, we used the form- 
factor model of Bergstrom et al. [l|, although we found 
little difference between results using it and a uniform 
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FIG. 5: The reconstructed pe + e~ mass distribution for the 
normalization mode after all cuts. The histogram is the sum 
of MC samples of Ep ee7 , K^ eeri decays and a uniform back- 
ground, where the relative amounts of each were determined 
by a fit, and the number of MC events was normalized to 
match the number of data events. The hatched area shows 
the main background source (uniform background). 
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TABLE I: Acceptances (A) and efficiencies (e) for the signal 
and normalization modes. 



Mode 


A(%) 


e{%) 


E+ MM decays 


0.259 


71.2 


Kp^ deca y s 


0.731 


69.4 


Normalization Ej ee7 


0.255 


5.6 



phase-space decay model. The form-factor model uses 
the SM processes (FCNC and internal conversion) shown 
in Fig. nja)-(c). In the model B{H + — > pl + l~) depends 
on three parameters: (62/61), (ci/61) and (C2/61). Us- 
ing the Particle Data Group values of the partial decay 
width and the a 7 decay parameter of the S + — > pry de- 
cay, we determined (62/61) = —0.46 ± 0.07. The other 
two parameters were chosen to make the £^~ ee branching 
ratio as small as possible so that it is consistent with the 
experimental upper limit, < 7 x 10~ 6 0.0]- This gives 
(ci/61) = 0.7 and (C2/61) = —13.0, and a branching ratio 
B(Y,+ -> pe+e~) = 8.9 x 10~ 6 . 

In the simulation of the SJp decay, we assumed 
that the P° mass was 214.3 MeV/c 2 with negligible de- 
cay width, and that the P° decayed immediately to the 
dimuon pair. The matrix element for the P° — > 
decay was taken to be uniform. 

MC simulations were used to estimate the geometric 
acceptances and the event-selection efficiencies for the 
signal and normalization modes. Table [i] shows a sum- 
mary of the results. The relative trigger efficiency for 
the signal mode with respect to the normalization mode 
and the unlike-sign dimuon detection efficiency were de- 
termined from the full data sample, and were found to 
be (91.9 ± 1.4)% and (96.2 ± 0.3)%, respectively. 

The systematic errors in the measurement of the 
branching ratios are listed in Table [H] They were stud- 
ied by varying the ranges of the kinematic parameters 
in the modeling of the S + production at the target, the 
assumed target positions, the magnetic fields, and the 
parameters describing the 5]+^ and ir® e ^ decays in the 
MC simulations. The systematic errors were dominated 
by uncertainties in the modeling of the S + production, 
which was mainly due to the small normalization sample. 
The systematic error includes the statistical error from 
the normalization mode. 

The branching ratio for the E^La decay hypothesis is 
[8.6ig;4(stat) ± 5.5(syst)] x 10~ 8 , in apparent disagree- 
ment with the expected branching ratio using the form- 
factor parameters given above, (1.1 ± 0.3) x 10~ 8 , where 
the error includes only the uncertainty in the parameter 
(62 /61 ) . The branching ratio for the ^pp MM hypothesis is 
£(£+ -^ P P°,P° -> (1+/J,-) = [3.1±^(stat)±1.5(syst)]x 
10~ 8 . The statistical errors for the branching ratio mea- 
surements were estimated by using the statistical table 
in Ref. |8j . If the three signal events are assumed to be 
background events from some unknown source, then us- 
ing the method in Ref. |9J to take the systematic error 



TABLE II: Fractional systematic errors (<jb/B) in the branch- 
ing ratios of E+ MM and of E+ Pflfl . 



Source 


OB 

y+ 


IB (%) 






E+ 


Normalization 


14.7 


14.7 


Modeling of E + production 


54.3 


44.6 


Beam targeting 


11.1 


8.7 


Magnetic field 


2.2 


3.9 


-Liiggei enioieiicy 


1.0 


1 . 


Muon identification 


0.3 


0.3 


E+j M form factor 


28.9 




7T° form factor 


1.8 


1.8 


B(E+ -> pit ) 


0.6 


0.6 


B(ty° -» e + e"7) 


2.7 


2.7 


MC statistics 


1.3 


1.3 


Total 


64.4 


48.1 



into account, we obtain an upper limit at 90% C.L. of 
B(Z+^) < 3.4 x 10- 7 . 

In summary, we observe three clean E + — > p\x + \i~ can- 
didates. This is the first evidence for this decay. The 
probability that the three events have a dimuon mass 
within 1 MeV/c 2 of each other in the SM is estimated 
to be < 1%. The three events are consistent with the 
E + — > pP°,P a — > decay hypothesis, with a P° 

mass of 214.3 ± 0.5 MeV/c 2 . 
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